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ABSTRACT 


Oblique  ionograms — records  of  sweep- frequency  ionospheric  sounders 
with  transmitter- receiver  separation  of  hundreds  or  thousands  of 
kilometers — differ  from  vertical-incidence  ionograms  in  that  they  cannot 
be  inverted  to  give  ionospheric  electron  density  without  making  unreal¬ 
istically  drastic  assumptions  concerning  the  nature  of  the  ionosphere. 
Oblique  ionograms,  however,  can  be  calculated  from  a  known  ionosphere. 
This  report  describes  a  method  for  performing  the  calculation  using  a 
digital  computer. 

The  technique  takes  full  account  of  earth  and  ionosphere  curvatures 
and  allows  arbitrary  choice  of  electron  density  with  height  and  longitu¬ 
dinal  distance  but  not  with  lateral  distance.  The  geomagnetic  field  is 
neglected  in  the  interest  of  economy,  but  it  could  be  included.  Examples 
of  synthetic  oblique  ionograms  for  a  Chapman  ionosphere  are  presented, 
both  with  and  without  an  E  layer,  over  path  lengths  from  0  to  8000  km. 
Results  included  show  the  effects  of  layer  height,  thickness,  tilts,  and 
a  variety  of  non-Chapman  electron-density  distributions.  The  synthetic 
ionograms  include  the  raypath  takeoff  angles  and  also  the  absolute  group 
time  delay;  these  are  particularly  useful  since  they  are  not  available 
on  experimental  ionograms. 
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I.  INTRODUCTION 


Many  research  organizations  throughout  the  world  are  presently  col¬ 
lecting  data  in  the  form  of  oblique  ionograms.  Such  data  are  directly 
useful  in  determining  the  communication  channels  that  are  open  along  the 
particular  path  spanned  by  the  oblique  sounders.  These  data  also  show 
great  promise  for  use  in  the  study  of  upper-atmospheric  electron-density 
distribution,  which  may  lead  to  a  better  understanding  of  the  behavior 
of  the  ionosphere.  The  primary  difficulty  at  the  moment  is  that  there 
is  no  practical  way  to  determine  the  electron-density  distribution  from 
an  oblique  ionogram  unless  one  is  willing  to  make  gross  simplifying 
assumpt ions . 

The  purpose  of  this  paper  is  to  introduce  a  technique  that  may  prove 
tc  be  of  significant  value  as  a  part  of  the  current  effort  to  gain  a 
better  understanding  of  oblique  ionograms.  If  one  begins  with  a  known 
a  priori  electron-density  distribution,  it  is  economically  feasible  to 
use  digital-computer  raytracing  techniques  for  the  synthetic  generation 
of  the  corresponding  oblique  ionograms.  The  basic  idea,  then,  is  that 
one  can  generate  a  library  of  synthetic  ionograms  whose  structure  and 
detail  can  be  analyzed  and  understood  completely.  Later  one  could  refer 
to  this  library  and  compare  the  synthetic  ionograms  with  actual  expert 
mental  ionograms  with  a  hope  that  the  comparison  will  aid  in  drawing 
uselul  conclusions  concerning  the  nature  of  electron-density  distribution 
which  actually  existed  along  the  path. 

The  synthetic  records  also  provide  information  not  available  on 
their  real  counterparts.  During  the  manufacture  of  the  synthetic  iono¬ 
grams,  it  is  particularly  easy  to  determine  the  takeoff  angles  of  the 
raypaths  that  lead  to  various  details  of  the  structure  of  the  oblique 
ionograms.  This  information  is  not  usually  available  and  it  should  be 
useful  to  those  who  design  antennas  used  in  gathering  oblique  ionograms. 
The  takeoff  angles  should  also  prove  helpful  in  the  interpretation  of 
actual  ionograms  because  they  show  which  portions  of  the  total  echo 
structure  are  rendered  invisible  because  of  low  antenna  gain  at  the 
particular  vertical  angles  in  question. 
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The  synthetic  records  show  the  actual  time  delay  between  transmitter 
and  receiver.  This  information  is  helpful  because  experimental  records 
have  arbitrary  time  base-lines  as  a  result  of  an  inherent  inability  to 
determine  the  total  accumulated  time  delay  accurately.  In  addition, 
synthetic  records  could  show  the  signal  strength  of  each  component  of 
each  trace,  but  these  data  were  not  included  on  the  records  given  here. 
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II. 


RAYTRACING  AND  THE  RAYSET  INFORMAT ION- STORAGE  TECHNIQUE 


A  large  amount  of  the  work  that  must  be  done  to  make  synthetic  iono- 
gi-ams  consists  of  the  tracing  of  radio  waves  in  the  upper  atmosphere  by 
a  digital  computer.  A  number  of  organizations  currently  have  raytracing 
programs  in  operation.  The  particular  program  used  for  the  generation 
of  data  in  this  paper  was  described  in  a  previous  report  by  the  author 
[Ref.  l],  It  should  be  emphasized  that  any  digital  raytracing  program 
could  be  used  in  exactly  the  same  way  as  described  here,  provided  that 
it  accurately  calculates  the  range  and  group  time  delay  of  each  ray. 

The  raytracing  program  used  here  consisted  basically  of  the  successive 
solutions  of  Snell's  Law  as  a  ray  was  followed  through  a  horizontally 
stratified  ionosphere  wherein  each  stratum  was  about  1  km  thick.  Because 
the  geomagnetic  field  was  neglected,  the  ionogrnms  generated  in  this 
study  will  be  found  to  lack  an  extraordinary  trace. 

One  of  the  major  factors  that  must  be  considered  in  synthesizing 
oblique  ionogrnms  is  the  cost  of  raytracing.  Unfortunately,  this  cost 
is  a  difficult  factor  to  evaluate  because  It  can  vary  by  an  order  of 
magnitude,  depending  on  the  particular  situation  at  hand.  Nevertheless, 
the  following  rough  estimates  are  offered:  for  each  hop  of  each  ray  in 
an  ionosphere  consisting  of  300  strata,  each  1  km  thick,  the  cost  is 
^0.01.  If  horizontal  gradients  of  electron  density  (tilts)  are  present 
in  the  ionosphere,  then  the  cost  goes  up,  to  approximately  50.10.  Although 
no  data  are  available,  it  is  this  author's  understanding  that  the  in¬ 
clusion  of  the  magnetic  field  drives  the  cost  up  by  another  factor  of 
perhaps  10,  so  that  the  cost  could  be  as  much  as  $1.00  per  hop  per  ray. 

The  importance  of  cost  control  can  be  appreciated  when  one  considers 

5 

that  this  lalroratory  alone  traces  approximately  10  rays  per  year. 

The  building  block  of  this  method  is  the  "rayset"— the  name  given 
to  the  set  of  numbers  that  characterizes  the  raypath  information  acquired 
through  raytracing.  Each  such  set  of  numbers  describes  one  hop  of  one 
ray  and  is  stored  on  one  punched  card.  For  example,  if  a  ray  executes 
three  hops,  the  computer  will  make  three  punched  cards  containing  all 
the  key  information  that  characterizes  each  hop.  This  system  is  illus¬ 
trated  by  Fig.  1. 
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KIG.  1.  THE  NATURE  OF  "RAYSETS . " 


Consider  the  process  that  goes  on  inside  a  digital  computer  during 
raytracing.  Each  ray  is  calculated  individually.  The  computer  is  in¬ 
structed  to  go  to  a  designated  starting  point  and  consider  a  ray,  origi¬ 
nating  in  that  point,  propagating  at  a  given  initial  takeoff  angle  above 
the  surface  of  the  earth.  The  computer  then  follows  that  ray  in  circular 
coordinates  and,  using  appropriate  equations,  it  calculates  the  curved 
path  of  the  ray  as  it  traverses  a  medium  of  varying  index  of  refraction. 
Many  things  can  happen  to  the  ray.  It  can  be  refracted  by  the  ionosphere 
and  execute  normal  hops  as  it  docs  in  the  first  two  hops  of  Fig.  1.  It 
can  go  so  far  that  it  exceeds  a  specified  maximum  range  and  the  computer 
will  be  instructed  to  shut  off.  Similarly,  the  ray  might  exceed  a 
specified  maximum  height.  With  horizontal  gradients  present,  the  ray 
may  take  off  and  climb  high  into  the  ionosphere  only  to  come  down  on  top 
of  a  densely  ionized,  localized  area;  then  it  may  refract  back  upward. 

In  such  a  case  it  is  said  that  the  ray  executed  a  "perigee."  This  fre¬ 
quently  happens  when  raytracing  is  conducted  in  an  ionosphere  with  a 
E  layer  whose  critical  frequency  increases  with  increasing  range. 
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It  may  happen  that  the  ray  trapped  between  the  E  and  F  layers 
will  eventually  escape.  However  the  valley  in  the  N-h  profile  between 
the  E  and  F  layers  may  fill  up  with  electrons  as  the  range  increases 
further,  and  the  ray  may  find  that  it  cannot  continue  to  propagate  be¬ 
cause  the  index  of  refraction  is  imaginary  in  all  forward  directions. 

In  such  a  case  the  ray  is  terminated  by  the  computer,  although  physically 
the  energy  must  reflect  back  out  of  such  a  closed  trap. 

As  the  computer  follows  the  progress  of  each  ray,  it  momentarily 
stops  every  time  a  ray  encounters  the  earth,  a  perigee,  or  its  termination. 
At  that  time  the  computer  stores  on  magnetic  tape  the  set  of  numbers, 
called  a  rayset,  that  will  later  be  punched  into  a  card. 

Consider  the  first  hop;  when  the  computer  finds  that  the  ray  strikes 
the  earth,  as  it  does  in  the  center  of  Fig.  1,  it  stores  the  numbers 
shown  in  the  top  row.  The  IID  (ionosphere  Identification)  is  a  three- 
digit  number  that  identifies  the  particular  electron-density  distribution; 
the  frequency  is  the  number  of  kilocycles;  the  parameter  B  (bounce) 
will  be  given  as  1  because  this  is  the  first  time  the  ray  has  hit  the 
earth;  the  parameter  0  is  the  takeoff  angle;  and  is  the  landing 

angle  as  shown.  The  three  parameters  Tgl>  Tr  and  Rj  are  the  group 

time  delay,  phase  time  delay,  and  the  ground  range  from  the  transmitter 
to  the  point  of  impact.  The  height  is  the  height  of  the  apogee  of 

the  ray  as  shown.  The  parameter  1^  (indicator)  is  a  number  that  is 
coded  0  through  6  and  tells  just  what  the  ray  did.  N'otice  that  I  =  0 
when  the  ray  hits  the  earth,  1=1  if  the  ray  exceeds  maximum  height, 
and  1=2  or  3  if  the  ray  exceeds  maximum  range.  Similarly,  I  =  A 
or  5  if  the  ray  executes  perigees  and  1=6  if  the  ray  encounters  a 
closed  trap. 

When  the  computer  detects  that  the  ray  has  encountered  the  ground 
lor  a  second  time,  ns  shown  at  the  right  of  lig.  1,  it  stores  the  second 
rayset  shown  across  the  center  of  the  figure.  Here  the  information  is 
similar  to  that  in  the  first  rayset  except  that  B  =  2  because  this  is 
the  second  encounter  with  the  ground. 

As  a  matter  of  interest,  it  might  be  noted  that  the  definitions  of 
the  various  parameters  as  shown  in  Fig.  1  apply  only  when  1=0.  When 
I  is  nonzero,  the  ranges,  times,  angles,  and  heights  are  defined 
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differently,  in  such  a  way  as  to  maximize  the  ability  of  the  user  to  ri 
construct  the  ray  trajectory.  As  an  example  of  this,  consider  that  the 
third  hop  of  a  ray  exceeds  maximum  range,  going  upward  as  it  appears  to 
do  on  Fig.  1;  the  indicator  would  be  1=2.  In  this  case,  the  time 
and  the  range  would  be  measured  to  the  point  of  termination,  the  height 
would  be  the  height  of  termination,  and  the  angle  \|r  would  be  the  angle 
of  the  ray  relative  to  the  horizontal.  If  the  ray  had  exceeded  maximum 
range  coming  down,  these  same  definitions  would  have  applied  except  that 
the  indicator  would  be  1=3. 
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III.  METHOD  OF  OBLIQUE- I ONOGRAM  SYNTHESIS 


In  order  to  make  an  artificial  oblique  ionogram  similar  to  one  that 
might  be  obtained  between  two  stations  3000  km  apart,  a  common  first 
reaction  is  an  inclination  to  compute  rays  of  several  frequencies  that 
propagate  a  distance  very  close  to  3000  km.  This  technique  would  yield 
the  desired  answers  and,  in  fact,  the  author  has  developed  a  specialized 
program  that  will  search  for  and  calculate  only  those  rays  that  go  to 
specified  ranges.  However,  it  turns  out  to  be  necessary  to  calculate 
an  entire  family  of  rays  during  the  search  process  when  looking  for  the 
takeoff  angle  that  will  lead  a  ray  to  the  desired  range.  Consequently, 
it  is  a  vcr\  inefficient  process  if  one  uses  information  onl y  from  rays 
that  propagate  the  desired  distance. 

The  method  that  has  been  used  to  produce  the  ionograms  given  here 
is  probably  near  to  being  the  most  economical  approach.  In  essence  the 
raysets  are  used  as  datn  for  a  digital  program  that  first  separates  the 
rayset  cards  so  that  it  may  consider  at  any  one  time  only  those  cards 
that  have  the  same  hop  number,  U.  Considering  one  such  group  of  cards 
and  one  frequency  at  a  time,  the  program  interpolates  in  a  three- 
dimensional  space  whose  axes  arc  ground  range  H,  group  t  ima  delay  T^, 
and  takeoff  angle  3.  The  program  considers  as  a  group  all  those  raysets 
that  have  the  same  B  and  frequency  F  and,  within  these  groups,  it 
considers  the  raysets  in  order  of  increasing  3.  In  this  order,  it  ex¬ 
amines  successive  values  of  K  until  it  finds  two  values  that  are  the 
end  points  of  an  interval  containing  the  desired  range.  It  then  assumes 
that  between  the  two  raysets  all  relations  arc  linear  and  it  performs 
an  interpolation  to  find  the  values  of  T  r  and  3  at  the  desired  range. 

In  the  majority  of  circumstances,  this  procedure  produces  answers 
that  arc  quite  reasonable,  since  for  small  increments  of  takeoff  angle 
the  relationships  really  arc  nearly  linear  between  successive  raypaths. 
On  occasion,  however,  this  assumption  breaks  down  and  it  is  necessary 
for  human  intervention  at  this  point-- examinat ion  of  the  printed  raysets 
and  comparison  with  the  computed  answers — to  sec  that  the  assumption 
is  truly  reasonable.  One  can  tell  from  trends  exhibited  by  the  numbers 
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In  cases  where  it 


whether  or  not  linearity  is  a  valid  approximation, 
appears  that  nonlinear  relationships  actually  must  hold,  it  is  neeessai> 
to  interpolate  according  to  some  logical  nonlinear  scheme,  or  else  to 
perform  more  ray  tracing. 

This  process  is  usually  carried  out  at  frequency  intervals  of 
approximately  4  Me.  The  answers  are  then  plotted  on  a  graph  of  group 
time  delay  vs  frequency  and  this  plot  _is  the  synthetic  ionogram.  After 
such  a  preliminary  plot  is  made,  the  one-hop  data  points  can  be  connected 
with  a  line,  and  the  two-  and  three-hop  data  can  be  treated  similarly. 

The  computer  tells  whether  each  point  is  an  upper  or  lower  ray,  as 
determined  from  the  algebraic  sign  of  dR/d3 .  From  these  lines  and  the 
blank  spaces  between  them,  the  exact  frequencies  and  takeoff  angles 
that  still  need  to  be  rayt raced  can  be  deduced  with  fairly  good  accuracy. 

Since  this  method  involves  interpolation,  the  synthetic  ionograms 
show  the  data  points  that  were  actually  calculated.  Thus,  the  reader 
can  judge  for  himself  which  of  the  results  presented  here  are  based  on 
mathematics  and  which  are  based  on  judgment. 

In  general,  if  the  electron-density  distribution  is  n  smooth  curve 
when  plotted  ns  a  function  of  height  and  range,  then  the  ionogram  will 
also  consist  of  smooth  curves.  In  this  situation,  n  synthetic  ionogram 
can  t>e  produced  with  a  minimum  of  raytracing  because  it  is  possible  to 
interpolate  over  greater  ranges  when  parameters  are  smoothly  varying 
functions.  Thus  the  amount  of  digital  calculation  can  be  kept  to  a  min¬ 
imum.  It  is  a  feature  of  the  method  of  simulation  that,  once  a  synthetic 
ionogram  is  made  for  a  given  ionosphere  and  a  given  range,  it  is  con¬ 
siderably  less  expensive  to  make  other  ionograms  for  different  ranges 
for  the  same  ionosphere.  This  advantage  results  because  the  raytracing, 
once  done,  can  be  used  again  in  the  inexpensive  "obi ique- ionogram  pro¬ 
gram"  at  new,  different  ranges.  The  only  cost  of  producing  such  extra 
ionograms  is  the  cost  of  the  few  special  rays  that  must  be  computed,  to¬ 
gether  with  the  cost  of  personnel  time. 
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IV. 


VERT  TPAT.-  INCIDENCE  IONOGRAMS 


A  vertical-incidence  ionogram  is  essentially  an  oblique  ionogram 
at  zero  range.  The  synthesis  of  these  records  is  extremely  easy  since 


it  is  necessary  only  to  integrate  group  time  delay  verti 
computer  program  was  written  to  carry  out  this  procedure 
same  ionosphere  description  cards  that  were  made  for  the 
program;  six  of  the  resulting  records  are  included. 


rally.  A  special 
,  utilizing  the 
raytracing 
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V.  SYNTHETIC  IONOGRAMS 


In  this  section  39  synthetic  oblique  ionograms  and  six  vertical 
ionograms  are  presented  that  were  generated  by  the  method  described. 
Since  this  is  the  first  effort  along  these  lines,  an  attempt  has  been 
mode  to  select  examples  of  electron-density  distribution  that  illustrate 
the  most  commonly  encountered  gross  characteristics  of  the  ionosphere. 

Nine  ionograms  are  made  with  the  Chapman  ionosphere  and  various 
ranges  between  transmitter  and  receiver.  Then  an  E  layer  is  added  to 
the  Chapman  layer  and  the  ionograms  are  recalculated.  The  effect  of  the 
density  of  the  E  layer  is  investigated  briefly.  Another  group  of  iono¬ 
grams  shows  what  happens  when  Chapman- layer  height  and  thickness  are 
changed.  Positive  and  negative  horizontal  gradients  of  electron  density 
are  studied  in  a  group  of  six  ionograms.  Then,  two  different  parabolic 
layers  are  investigated  at  each  of  five  ranges.  Finally,  a  rather  un¬ 
usual,  layered  ionosphere  1  used  to  generate  some  fairly  intricate 
ionograms . 

When  a  Chapman  layer  is  mentioned  in  the  following  descriptions, 
this  notation  is  implied: 

N  sr  electron  density  fin  electrons/cc) 

\'0  =  maximum  electron  density  in  the  layer 

11  =  height  of  maximum  electron  density  N  (in  km) 

m  u 

H  =  scale  height  of  Chapman  layer  (in  km) 
s 

7.  =  (hcight-llm)/(i!s) . 

With  this  notation,  a  Chapman  layer  means  an  electron  density  N 
such  that 


S  =  Nq  exp( 1  -  z  -  c  /)  . 

Before  the  data  are  presented,  one  further  item  deserves  mention. 
Because  this  work  neglects  the  magnetic  field,  tnc  actual  magnitude  of 
the  frequencies  at  which  calculations  are  made  is  not  really  the  variable 
of  interest.  The  more  significant  variable  is  the  ratio  of  the  radio 
frequency  to  the  plasma  frequency  of  the  medium.  This  conclusion  will 
be  seen  to  follow  if  one  considers  that  the  index  of  refraction  used  in 
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computations  of  ray  paths  is  simply  n, 


where 


n2  =  1  -  (f 2Ji2)  , 

in  which  f  is  the  local  plasma  frequency  and  f  is  the  operating 
frequency.  Therefore,  if  following  results  indicate  that  a  certain  time 
delay  exists  along  a  certain  path  at  some  specific  frequency,  then  the 
time  delay  presented  would  also  be  correct  for  half  that  frequency, 
provided  the  ionospheric  electron  density  were  reduced  throughout  the 
entire  range  of  heights  by  the  factor  (l/2)2.  Any  other  factor  will 
also  work. 

As  n  consequence  of  this  reasoning,  it  might  be  convenient  to  plot 
both  the  oblique- incidence  ionograms  and  their  associated  electron- 
density  profiles  on  graphs  such  that  the  electron-density  and  frequency 
scales  are  logarithmic.  Then  it  would  not  be  necessary  to  commit  oneself 
to  a  specific  frequency  scale  on  the  ionograms  but,  rather,  one  could 
label  the  horizontal  axis  in  terms  of  the  parameter  f/f^  where  fp 
is  the  critical  frequency. 

Most  experimental  ionograms,  however,  arc  presently  being  made  by 
the  equipment  manufactured  by  Granger  Associates,  wherein  the  frequency 
scale  is  a  sequence  of  four  different  linenr  scales.  Ir.  an  effort  to 
maximize  the  similarity  of  these  synthetic  ionograms  to  the  actual  experi¬ 
mental  data  that  are  available,  they  arc  presented  on  a  "Granger  Scale," 
which  makes  it  necessary  to  assign  specific  frequencies  to  the  horizontal 
axis.  Each  such  plot  will  be  accompanied  by  a  corresponding  ionogram 
that  has  a  logarithmic  frequency  scale,  since  some  oblique  soundcis 
present  their  data  in  this  manner.  Also,  this  approach  will  permit  the 
reader  to  apply  the  si iding-frequcncy-scalc  principle  if  he  so  desires. 

Small  solid  triangles  along  the  ionogram  frequency  axes  indicate 
the  vert ical-inciUence  critical  frequency  of  the  particular  ionosphere. 

In  most  cases  this  value  is  near  9  Me,  because  most  of  the  synthetic 
ionospheres  have  arbitrarily  been  assigned  to  a  maximum  electron  density 
of  10G/cc.  Also,  the  small  open  triangles  in  the  body  of  the  ionogram 
indicate  frequencies  at  which  raytracing  was  done  without  finding  any 
propagation  modes  at  the  particular  range. 
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The  format  of  the  vert ical- incidence  ionogrnms  was  chosen  to  maxi¬ 
mize  its  similarity  to  the  standard  format  with  which  this  author  is 
familiar.  Also,  the  vert ical- ionogram  records  show  the  t rue-height 
profile,"  which  is  simply  the  input  data  describing  the  vertical  distri¬ 
bution  of  electron  density.  This  profile  is  interesting  in  that  it  shows 
the  answer  which  one  should  get  from  a  reduction  of  the  vertical- 
incidence  record. 


A.  A  CHAPMAN  LAYER  AND  VARIOUS  RANGES 

Figure  2  shows  the  electron-density  distribution  in  the  Chapman 
layer  whose  I  ID  ( Ionosphere  Ident if lent  ion)  is  Oil .  This  ionosphere 
has  „m  =  300  km,  l!s  =  100  km,  and  NQ  =  106.  These  parameters  were 
chosenVeause  they  produce  a  density  that  is  typical  of  the  actual  dis¬ 
tribution  encountered  in  practice.  Later,  the  effect  of  vaiintions 
these  parameters  will  be  illustrated.  Notice  that  fc  =  9  Me,  approxi¬ 
mately,  but  that  this  value  need  not  be  regarded  ns  a  fixed  quantity  be¬ 
cause  of  our  ability  to  apply  a  sli-  ing-scalc  concept  to  the  plasma  and 
operating  frequencies,  as  previously  mentioned. 

Figures  3  through  11  show  ionogrnms  that  were  synthetically  gen¬ 
erated  with  the  assumption  that  the  separation  distances  between  the 
transmitter  and  receiver  were  0,  1000.  2000,  3000,  1000,  5000,  6000, 

7000,  and  8000  km,  in  that  sequence.  Notice  that  under  this  scheme, 

Fig.  3  is  simply  a  vert ical- incidence  ionogram;  Fig.  1  is  an  oblique 
ionogram  in  which  the  stations  are  separated  by  1000  km,  and  so  forth. 
Each  trace  is  labeled  to  indicate  the  number  of  hops  of  the  rays  that 
generated  it.  The  upper-  and  lower-ray  portions  of  each  echo  are 
easily  identified  and  thus  are  not  labeled.  Each  point  that  has  been 
calculated  is  identified  by  a  dot  along  the  curve  on  the  "Granger"  plot. 
Near  most  of  the  dots  are  small  numbers,  which  indicate  the  takeoff 

angles  of  the  particular  rays  involved. 

For  these  and  for  most  of  the  ionogrnms  in  this  report,  calculation 
was  cut  off  at  three  hops;  therefore,  the  absence  of  four-hop  returns 
is  not  significant.  Because  this  particular  ionosphere  does  not  have 
any  horizontal  electron-density  gradients,  it  would  be  possible  to  gen¬ 
erate  a  fourth-hop  return  on  the  8000- km  plot  by  simply  doubling  the 
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numbers  associated 
such  a  process  can 
decision  was  taken 


with  the  two-hop  return  on  the  4000-km  plot.  However, 
be  carried  on  without  end  and  consequently  the  arbitrart 
to  stop  at  three  (or  later  four)  hops. 
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FIG.  2.  ELECTRON- DENSITY  DISTRIBUTION  FOR  CHAPMAN  LAYER,  IID  Oil. 
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FREQUENCY  (Me) 


0  I  I  at  1000  km 


FIG.  4.  I0N0GRAM  FOR  IID  Oil;  SEPARATION  DISTANCE,  1000  km. 
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FIG.  5.  IONOGRAM  FOR  IID  Oil;  SEPARATION  DISTANCE,  2000  km. 
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Oil  at  3000  km 


fig.  g.  ionogram  FOR  I id  oil;  separation  DISTANCE,  3000  km. 


1  R 


SEL- 64-106 


granger 

FREOUENCY 

SCALE 


POSSIBLE  TIME,  R/C 


7  8  10 


logarithmic 

frequency 

SCALE 


, _ , - - — i — i — , — i — i — i — i — i — i — iii  i  ii  .  r- 

T  20  24  28  32  40 


12  14  16 

FREOUENCY  (Me) 


Oil  at  4000  km 

IONOGRAM  FOR  I ID  Oil;  SEPARATION  DISTANCE,  4000  km. 
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FIG.  9. 


IONOGRAM  FOR 


I  ip  oil;  SEPARATION  DISTANCE,  6000  km 
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FIG.  11.  IONOGRAM  FOR  1ID  Oil;  SEPARATION  DISTANCE,  8000  km 
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B.  WEAK  AND  STRONG  E-LAYERS  AT  VARIOUS  RANGES 


In  order  to  investigate  the  effect  of  an  E  layer  on  oblique  lono- 
grams,  the  elctron-density  distribution  IID  Oil  was  modified  with  a 
weak  E  layer  ( IID  012)  and  a  strong  E  layer  ( IID  089)  as  shown  in 
Fig.  12.  These  particular  electron  distributions  were  chosen  subjectively 
with  the  idea  that  012  is  fairly  typical  and  089  is  more  prominent  than 

the  average  E  layer  encountered  in  nature. 

Ionograms  for  012  were  calculated  at  the  same  nine  ranges  presented 
on  the  Chapman  layer  Oil  (Figs.  13  through  21 ),  whereas,  because  089  is 
atypical,  its  ionogram  was  calculated  only  at  0  and  3000  km  (Figs.  22 
and  23).  Since  012  is  identical  to  011  in  all  respeets  except  the  ad¬ 
dition  of  the  E  layer,  the  comparisons  of  these  ionograms  are  partic¬ 
ularly  interesting.  One  can  also  compare  the  3000-km  ionograms  for  Oil, 
012,  and  089  and  see  the  effect  of  the  degree  of  ionization  in  the  E 

1 ayer. 

The  ionograms  calculated  for  Chapman  layer  Oil  are  repeated  in  red 
on  the  lower  (logarithmic)  presentations  or  these  ionograms  to  permit  a 
clearer  visualization  of  the  changes  caused  by  the  presence  of  the  E 
layer.  In  addition,  the  089  oblique  ionogram  shows  (in  blue)  the  cor¬ 
responding  trace  obtained  for  IID  012. 
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IONOGRAM  FOR  I  ID  012;  SEPARATION  DISTANCE, 
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012  at  1000  km 


IONOGRAM  FOR  IID  012;  SEPARATION  DISTANCE,  1000  km. 
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FIG.  15.  IONOGRAM  FOR  IID  012;  SEPARATION  DISTANCE,  2000  km. 
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012  at  5000  km 


IONOGRAM  FOR  IID  012;  SEPARATION  DISTANCE,  5000  km. 
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FIG.  20.  IONOGRAM  FOR  IID  012;  SEPARATION  DISTANCE,  7000  km 
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012  at  8000  km 

FIG.  21.  IONOGRAM  FOR  I  ID  012;  SEPARATION  DISTANCE,  8000  km. 
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C.  CHANGES  IN  LAYER  HEIGHT 

The  effect  of  changes  in  the  height  of  the  maximum  electron  density 
in  a  Chapman  layer  is  examined  here.  The  oblique  ionograms  are  calculated 
only  for  a  separation  distance  of  3000  km.  Two  new  ionospheres  were  gen¬ 
erated  with  H  's  selected  so  that  their  heights  bracketed  the  value  used 
’  m 

for  I  ID  Oil.  Also,  in  order  to  maintain  the  similarity  to  Oil,  it  was 
decided  to  maintain  the  ratio  \lJHs  =  3.  This  ratio  resulted  in  the 
following: 


IIP 

H  (  km) 
m 

Hs(km) 

053 

200 

66.7 

011 

300 

100 

054 

400 

133.3 

The  maximum  electron  density  for  the  new  ionospheres  is  still  10  per  cc, 
as  it  was  for  Oil.  These  new  electron-density  distributions  are  shown 
together  with  Oil  in  Fig.  24 .  The  3000-km  ionograms  calculated  using 
these  new  ionospheres  arc  shown  on  Figs.  23  and  26.  For  comparison  pur¬ 
poses,  the  appropriate  oblique  ionogram  for  ionosphere  Oil  is  repeated 
in  red  on  the  logarithmic  presentations  that  follow. 
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fig.  26.  IONOGRAM  FOR 
HIGH  CHAPMAN  LAYER. 


I  ID  054;  SEPARATION  DISTANCE,  3000  km; 
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D.  CHANGES  IN  SCALE  HEIGHT 

The  effect  of  changes  in  the  scale  height  of  a  Chapman  layer  (when 
the  height  of  the  maximum  density  is  fixed)  is  now  examined.  In  order  to 
do  this,  two  new  ionospheres  were  generated,  again  bracketing  the  stan 
dard"  Chapman  layer.  The  values  chosen  were  as  follows: 


IID 

H  (  km) 
m 

II  (km) 

s 

055 

300 

80 

Oil 

300 

100 

056 

300 

120 

These  ionospheres  are  shown  on  Fig.  27  for  comparison.  The  new  iono- 
grams,  generated  for  a  3000-km  separation  distance  only,  are  shown  in 
Figs.  28  and  29.  Again,  these  figures  show  the  corresponding  oblique 
ionogram  for  IID  Oil  in  red  on  the  logarithmic  presentations. 


41 


SEL-64-106 


E.  HORIZONTAL  GRADIENTS  WITHIN  A  CHAPMAN  LAYER 


To  investigate  the  effect  of  horizontal  gradients  on  synthetic 
ionograms,  a  tilted  ionosphere  was  generated  according  to  the  following 
specifications:  directly  above  the  transmitter,  the  electron-density 

distribution  varies  with  altitude  like  that  of  IID  Oil,  except  that 
N  =  2  \  106,  exactly  twice  the  electron  density  that  was  present  in 
HD  Oil.  If  one  traveled  4000  km  from  the  transmitter  in  the  direction 
of  the  receiver  and  then  measured  the  electron-density  distribution,  one 
would  again  find  a  Chapman  layer,  except  that  then  NQ  =  A  10  Thub 
the  electron-density  distribution  is  fixed  above  two  points  along  the 

path . 

In  order  to  specify  the  method  of  variation  from  one  extreme  to  the 
other,  it  was  decided  to  make  electron  density  a  linear  function  of 
range  between  0  and  4000  km  and  then  nontilted  beyond  4000.  This  rather 
complicated  distribution  is  labeled  IID  =  093.  A  similar  ionosphere 
(IID  094)  was  generated,  but  in  it  the  electron-density  profiles  at  the  two 
end  points  of  the  path  were  switched.  Another  way  of  looking  at  this  is 
to  think  of  a  1000- km  separation  distance  between  transmitter  and  receiver 
with  the  tilted  ionosphere  existing  only  between  them.  Then,  for  IID 
093,  the  transmitter  is  under  the  most  dense  ionosphere  and  the  receiver 
is  under  the  least  dense  ionosphere.  For  IIP  091,  the  transmitter  and 
receiver  positions  are  switched. 

1,  is  interesting  to  note  that  reciprocity  should  apply,  and  indeed 
il  does,  lor  if  the  reader  will  examine  the  1000-km  ionograms,  he  will 
see  that  they  are  indeed  identical,  even  though  different  ray-tracings 
were  used  to  generate  them.  The  takeoff  angles,  however,  are  different 
because  the  takeoff  angles  Tor  one  ionogram  are  the  landing  angles  of 

the  other. 

The  electron-density  distributions  of  the  end  points  of  these  two 
ionograms  arc  shown  on  Fig.  30.  Vert ica 1- incidence  ionograms  are  not 
shown  because  they  would  be  identical  in  form  to  those  of  I-ig.  3,  except 
for  this  scale  shift,  which  has  been  discussed.  For  each  of  these  two 
ionospheres  three  ionograms  arc  given,  in  Figs.  31  through  36. 
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The  logarithmic  presentation  of  each  ionogram  also  shows,  in  red, 
the  ionogram  that  would  be  obtained  if  one  calculated  on  the  basis  of 
the  electron-density  profile  above  the  path  midpoint  without  any  horizon¬ 
tal  tilt. 
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FIG.  30.  ELECTRON  DENSITIES  AT  PATH  ENDPOINTS  FOR  IIDs  093  AND  094. 
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FIG.  31.  IONOGRAM  FOR  I ID  093;  SEPARATION  DISTANCE,  2000  km. 
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HG.  33.  IONOGRAM  FOR  I ID  093;  SEPARATION  DISTANCE,  6000  km 
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FIG.  34.  I0N0GRAM  FOR  IID  094;  SEPARATION  DISTANCE,  2000  km. 
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094  at  4000  km 


FIG.  35.  IONOGRAM  FOR  I  ID  094;  SEPARATION  DISTANCE,  4000  km. 
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F.  TWO  PARABOLIC  LAYERS  AT  VARIOUS  RANGES 

Parabolic  layers  are  of  interest  in  raytracing  primarily  because 
they  constitute  a  special  case  in  which  raypaths  can  be  solved  analyt¬ 
ically  without  the  use  of  a  digital  computer  or  its  associated  limitations, 
Consequently,  a  parabolic  layer  can  be  put  into  the  digital-computer 
program  and  the  results  can  be  checked  by  analytic  methods.  However,  it 
should  be  appreciated  that  in  nature  the  electron-density  distribution 
is  not  very  similar  to  any  parabolic  layer  except  perhaps  near  the  nose 
of  the  curve;  the  usefulness  of  the  parabolic  concept  is  primarily  math¬ 
ematical  . 

Some  time  ago  the  author  was  requested  to  make  synthetic  ionograms 
for  two  particular  parabolic  layers  at  four  specified  ranges.  The  results 
are  presented  here  because  they  are  of  fairly  general  interest  and  be¬ 
cause  they  do  show  what  oblique  ionograms  would  look  like  if  parabolic 
layers  actually  did  exist.  Also,  others  may  wish  to  use  them  to  check 

analytical  techniques  for  raytracing. 

The  two  parabolic  layers,  IID  111  and  12-1,  are  shown  on  Fig.  37, 
together  with  the  various  parameters  that  characterize  them.  Each  has 
been  used  to  calculate  synthetic  ionograms  nt  rnngcs  of  0,  1375,  1833, 

2750,  and  5500  km  to  four  hops.  These  results  arc  presented  on  ligs. 

38  through  47.  The  ionograms  cnlculntcd  for  IID  111  are  repeated  in 
red  on  the  logarithmic  presentations  of  the  1 2-1  oblique  ionograms  in 
order  to  facilitate  comparison.  It  can  be  seen  that  the  differences  be¬ 
tween  these  ionograms  are  due  solely  to  an  increase  in  the  height  of  the 
124  layer. 
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G.  37.  ELECTRON  DENSITIES  IN  TWO  PARABOLIC  LAYERS— IIDs  111 
AND  124. 
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FIG.  38.  IONOGRAM  FOR  IID  111;  SEPARATION  DISTANCE, 
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FIG.  39.  IONOGRAM  FOR  I  ID  111;  SEPARATION  DISTANCE,  1375  km. 


57 


SEL- 64-106 


fig.  40.  IONOGRAM  FOR  IID  111;  SEPARATION  DISTANCE,  1833  km. 
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HG.  41.  IONOGRAM  FOR  IID  111;  SEPARATION  DISTANCE,  2750  km. 
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III  at  5500  km 

FIG.  42.  I0N0GRAM  FOR  IID  111;  SEPARATION  DISTANCE,  5500  km. 
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FIG.  44.  IONOGRAM  FOR  IID  124;  SEPARATION  DISTANCE,  1375  km. 


62  - 


SEL- 64-106 


IONOGRAM  FOR  I ID  124;  SEPARATION  DISTANCE,  27dO  km 
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FIG.  47.  IONOGRAM  FOR  IID  124;  SEPARATION  DISTANCE,  5500  km. 
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G.  AN  UNUSUAL  EXAMPLE 

At  the  time  the  ionograms  just  presented  were  requested,  a  very 
unusual  example  was  also  requested;  it  is  presented  here.  Although 
the  electron-density  distribution  is  irregular,  the  features  found  on 
the  resulting  ionograms  are  exaggerated  examples  of  features  that  seem 
to  be  present  on  actual  oblique  ionograms. 

The  electron-density  distribution  shown  on  Fig.  48  was  generated 
from  sections  of  parabolic  layers.  The  resulting  five  ionograms  at 
various  ranges  are  presented  in  Figs.  49  through  53.  In  particular, 
notice  the  behavior  of  the  lower  raytraces  in  the  upper  layer.  The  ex¬ 
aggerated  way  in  which  the  range  of  these  traces  decreases  ns  the  fre¬ 
quency  increases  seems  to  be  characteristic  of  some  echoes  seen  on  ex¬ 
perimental  traces  and  is  apparently  due  to  the  long  interval  of  height 
in  which  there  is  no  change  in  electron  density.  The  phenomenon  is  due 
to  Pederson  rays  in  the  lower  layer  which  arc  reflected  from  the  upper 
layer.  As  the  frequency  increases,  the  retardntion  decreases  but  the 
ray  trajectories  remain  relatively  stable  with  the  consequent  apparent 
decrease  in  the  range  of  the  echoes. 
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FIG.  19.  IONOGRAM  FOR  IID  133;  SEPARATION  DISTANCE, 
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IONOGRAM  FOR  IID  133;  SEPARATION  DISTANCE,  1375  km. 
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FIG.  51.  IONOGRAM  FOR  I ID  133;  SEPARATION  DISTANCE,  1833  km. 
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VI. 


examples  showing  details  of  synthesis 


To  provide  some  insist  into  the  reason  for  various  features  on  the 
ionograms,  the  three  brief  sections  of  this  chapter  show  the  nature  of 
the  calculation  that  resulted  in  some  of  the  structures  presented.  The 
first  example  is  classically  simple,  the  second  example  is  more  complex 
_ _ _  and  the  third  points  out  an  unexpected  phenomenon. 


A.  DETAILS  OF  THE  NOSE  OF  I  ID  111  TWO-HOP  ECHO 


Reference  to  Fig.  <10  will  show  that  the  two-hop  echo  has  a  very  flat 
nose.  The  small  open  triangles  at  11.6  and  11.75  Me  show  that  raytracing 
was  done  and  that  no  propagation  path  existed  at  those  two  frequencies. 
The  situation  is  more  clearly  revealed  by  Fig.  54,  in  which  the  takeoff 
angles  and  two-hop  ground  ranges  of  the  various  pertinent  rays  are  shown. 
Notice  the  vertical  line  through  a  ground  range  of  1833  km.  Propagation 
paths  are  known  to  exist  at  that  range  whenever  the  curves  on  the  graph 

intersect  the  vertical  line. 
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The  search  process  during  the  generation  of  Fig.  40  involved  the 
following:  first,  an  11-Mc  raytracing  was  run,  which  revealed  the  ex¬ 

istence  of  a  propagation  path  at  a  takeoff  angle  near  39  deg.  From  this 
point  and  the  other  points  that  had  so  far  been  found  on  the  ionogram, 
it  could  be  determined  that  a  nose  of  the  two-hop  curve  should  be  nearby. 
Consequently,  raytracing  was  carried  out  at  11.50  and  11.75  Me.  As  seen 
in  Fig.  40,  the  11.50-Mc  raytracing  yielded  two  points  differing  in 
group  time  delay  by  a  little  over  3/4  msec  while  the  11.75-Mc  raytracing 
indicated  that  the  two-hop  skip  distance  was  near  1900  km.  Then,  from 
a  sketch  similar  to  Fig.  54,  it  was  determined  that  the  nose  of  the  curve 
should  be  very  near  11.6  Me.  The  raytrncing  was  done  as  shown  and, 
although  no  propagation  path  existed  to  1833  km,  the  lesults  were  suffi 
ciontly  close  to  enable  completion  of  the  nose  of  the  ionogram  with  a 
high  degree  of  confidence  in  the  resulting  shape. 


B.  DETAILS  OF  THK  E- LAYER  EFFECT  IN  II D  012 


This  section  will  demonstrate  in  some  depth  the  method  of  making 
the  synthetic  ionograms,  and  it  should  aid  in  understanding  the  reason 
lor  some  of  the  details  on  an  oblique  ionogram  when  E  layers  are  present. 
To  be  particularly  lucid,  a  small  group  of  frequencies  and  tnkeoff  angles 
was  selected  and  much  more  raytracing  was  done  with  this  group  than  was 
actually  necessary  to  produce  the  synthetic  records  presented  here.  This 
exercise  also  serves  to  show  that  the  output  of  the  computer  is  self- 
cons  1st  ent . 

Examination  of  the  ionograms  shows  that  there  is  a  great  deal  of 
structure  on  the  012  ionosphere  records  in  the  vicinity  or  15  Me.  Con¬ 
sequently,  raytracing  was  done  every  half  megacycle  between  14  and  17  Me. 
Each  ray  was  computed  to  three  hops,  resulting  in  357  raysets  distributed 
among  the  seven  frequencies,  and  cost  a  total  of  53.19.  The  51  raysets 
at  16.5  Me  are  printed  out  in  Fig.  55  to  illustrate  the  nature  of  these 
data.  (Costs  would  have  been  nearer  520  if  the  ionosphere  had  been 
tilted.) 

The  357  punched  cards  containing  raysets  were  then  processed  b\ 
another  computer  program,  which  carried  out  the  interpolation  required  to 
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FIG.  55.  PRINTED  RAYSETS  AT  16.5  Me  FOR  I ID  012. 
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give  all  the  obi ique- ionogr am  data  for  ranges  of  1000,  2000,  3000, 

5000  0000,  7000,  and  8000  to.  An  example  of  the  output  of  this  par,  c- 

ular' program  Is  shown  on  Fig.  56.  Each  lino  in  that  printout  represents 
a  propagation-path  delay  that  should  he  entered  on  one  o,  the  syn.hc 
ionograms.  For  example,  the  firs,  line  shows  that,  a,  frepuency  10.5 
on  HD  012,  there  Is  a  one-hop  upper-ray  mode  that  reaches  a  range 

,  a  r  i  an 9  rinir  and  has  an  excess  group  time 
3000  km,  has  a  takeoff  angle  of  1.802  deg,  ana  n 

delay  of  0.227  msec.  This  "excess  group  time  delay  (svmboli/cd 
1,  simply  the  accumulated  group  time  delay  along  the  3000-k,  path  fro. 
which  the  computer  has  subtracted  the  time  required  for  a  light  r»>  to 
traverse  3000  km.  It  turns  ou,  to  be  quite  convenient  to  work  .1, .  ex 

cess  delay  rather  than  total  delay,  because  the  excess  delay  is  nol  a 

„  The  last  column  gives  a  relative  measure 

strong  function  of  ground  range.  The  last  coium 

..  -  is  os  rib  along  this  particular  rnypath.  This 

of  the  signal  strength  ns  15.98  db  along 

information  was  not  used  on  ionograms  given  here. 

The  16. 5- Me  data  are  particularly  interesting  because  of  the  ray 

that  takes  off  at  2  deg.  This  ray  reaches  a  two-hop  range  of  6999.4652 

The  rays  both  above  and  below  this  particular  ray  reach  a  shorter 

.  tho  computer  gives  no  answer  for  7000  km, 
ground  range.  Consequently,  the  compute  g 

although  a  human  examiner  can  mahe  the  logical  Judgm.nl  that,  for  praclica 
purposes,  the  2-deg  ray  actually  did  reach  7000  k,  and  thus  should  be 
entered  on  the  oblique  lonogram.  This  Is  a  good  example  of  one  o  the 
many  reasons  lor  having  some  lot,  or  . . .erven,  lot,  in  long.  Into  «e 

sequences  of  computer  calculations. 

Seven  pages  of  calculations  such  as  those  shown  on  Fig.  *  we  e 

carried  out  by  an  IBM  7090  at  a  total  cost  of  11.25  (rental  i*»  -  ■> 

Figure  57  shows,  In  some-hat  idealized  form,  the  ranges  of  the  va. 
iouH  raypaths  that  were  computed  as  a  function  of  their  takeoff  angle 
,t  each  of  the  seven  frequencies.  Dotted  lines  are  shown  where  ex. c 

shape  of  the  crossover  curve  is  not  known.  The  reader  can  compare  1  g. 

55  with  the  16.5-Mc  curve  on  Fig.  57  to  see  how  the  lntt.r  wa  M  ^  . 

On  Fig .  55  it  can  be  seen  that  the  jump  in  range  between  1 . o-  and 
Leoff  angles  is,  indeed,  abrupt.  Figure  57  is  partitulatly  ntere  ting 

because  it  can  be  interpreted  to  show  the  basic  reason  for  most  of 
structures  on  the  HD  012  oblique  ionograms  at  frequencies  between 
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8000.0 
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FIG.  56.  COMPUTED  IONOGRAM 

POINTS 

TAKEN  FROM 

DATA  OF 

FIG.  55. 

and  17  Me.  For  example,  if  the 

6000-km 

ionogram  is  to  be 

generated. 

draw  a  vertical  line  through  3000  km;  the  intersections  of 

the  various 

curves  with  this  vertical  liie  would  indicate  the  existence  of  two-hop 
propagation  modes.  (This  works  only  lor  nontillod  ionospheres.)  Sim¬ 
ilarly,  one  could  draw  a  vertical  line  through  2000  km,  and  intersections 
with  this  line  would  represent  propagation  modes  via  three  hops  to  a 
ground  range  of  6000  km.  Thus,  if  one  draws  these  two  lines,  each  inter¬ 
section  will  represent  an  entry  of  data  on  the  ionogram.  This  process 
is  illustrated  in  Fig.  58,  where  the  6000-km  ionogram  is  shown  with  these 
intersections  plotted.  Notice  that  the  curves  are  smooth  except  for  the 
break  at  16  Me,  which  is  introduced  by  the  change  in  the  horizontal-scale 

compression  ral io  on  Granger  equipment. 

In  essence,  the  dots  on  Fig.  58  represent  the  amount  of  information 
that  a  digital  computer  can  easily  produce  without  human  intervention. 

Ii  price  were  no  object,  one  can  easily  see  that  the  entire  process  could 
be  made  totally  automatic. 
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ONE-HOP  GROUNO  RANGE  (km) 


FIG.  57.  GRAPH  OF  RAYPATH  RANGE  VS  TAKEOFF  ANGLE  TOR  IID  012. 
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KIG.  58.  PARTIALLY  SYNTHESIZED  IONOGRAM  SHOWING  POINTS  PLOTTED 
DIRECTLY  FROM  THE  COMPUTER  OUTPUT- 


40 
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c. 


ORIGIN  OF  A  ONE-HOP  TRACE  IN  A  5000-km  IONOGRAM 


\ 


A  commonly  encountered  "rule  of  thumb"  says  that  one-hop  rays  travel 
only  about  4000  km.  Consequently,  it  is  rather  surprising  to  find  a  one- 
hop  trace  on  a  5000-km  ionogram  such  as  is  shown  for  IID  Oil  on  Fig.  8. 
This  phenomenon  can  be  visualized  more  clearly  using  Fig.  59,  which 
shows  the  takeoff  angle  of  various  rays  plotted  against  the  ground  range 
that  the  rays  achieved  in  one  hop.  Notice  that  the  energy  density  at  the 
ground  will  be  proportional  to  the  quantity  d3/dR,  which  is  the  slope 
of  the  various  curved  lines  on  the  figure.  (However,  energy  density  does 
not  become  infinite  when  the  curves  arc  vertical;  see  for  example,  Ref.  2.) 

A  respectably  high  energy  density  propagates  to  ranges  well  in  excess 
of  5000  km.  In  fact,  it  can  be  seen  that  strong  31-Mc  signal  levels 
exist  between  5460  and  5500  km  via  lower  rays.  The  cause  of  the  one-hop 
trace  on  Fig.  8  can  be  understood  by  examination  of  the  intersections  of 
the  curved  lines  with  the  vertical  line  through  5000  km.  It  can  be  seen 
that  the  echo  has  significant  strength  only  in  the  frequency  interval  of 
approximately  30.8  to  30.9  Me,  and  these  signals  propagate  at  extremely 
low  takeoff  angles,  which  are  normally  heavily  attenuated  in  the  terminal 
equipment . 

Also,  it  should  be  pointed  out  that  this  electron-density  distribution 
has  a  far  smoother  N-h  profile  than  is  normally  encountered  in  nnture, 
and  the  resulting  families  of  mathematically  smooth  curves  such  as  on 
Fig.  59  probably  do  not  exist  in  the  natural  situation.  Experimental 
oblique  soundings  showing  one-hop  propagation  at  a  range  of  5600  km  arc 
shown  and  analyzed  in  Ref.  3. 
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KIG.  59.  I  XING- RANGE,  ONE-HOP  PROPAGATION  IN  A  CHAPMAN  1JVYEK  FOR  I  ID  Oil. 
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VII.  SELECTED  EXPERIMENTAL  IONOGRAMS  SHOWING  AGREEMENT  WITH  THEORY 


N* 


A  few  experimental  oblique  ionograms  are  presented  here  to  show 
that  the  structures  that  have  been  synthesized  are  actually  found  in 
nature.  All  of  the  records  given  here  were  selected  from  Refs.  4,  5,  and 
6,  and  are  reproduced  with  the  kind  permission  of  the  authors.  A  more 
thorough  literature  search  would  undoubtedly  turn  up  many  more  records 
exhibiting  appropriate  structures. 

On  Fig.  60  are  a  number  of  examples  illustrating  the  situation  when 
the  electron-density  profile  is  a  single  smooth  curve,  such  as  a  Chapman 
layer  (IID  Oil,  053,  054,  055,  056)  or  a  parabolic  layer  (HD  111  or 
124).  In  the  original  presentations  the  vertical  axes  were  labeled  in 
kilometers,  but  two  different  conversion  constants  were  used:  300  km  per 
msec,  and  150  km  per  msec.  Unfortunately  both  of  these  conventions  can 
be  supported  by  a  line  of  reasoning.  In  order  to  avoid  confusion,  all 
vertical  axes  here  are  labeled  in  msec,  which  is  unambiguous.  Also,  the 
figures  from  Ref.  4  have  a  logarithmic  frequency  scale  whereas  the  figures 
from  Refs.  5  and  6  have  a  linear  frequency  scale. 

Notice  the  striking  similarities  between  these  ionograms  and  those 
that  were  presented,  for  example,  in  the  first  set  of  ionograms  of  IID 
011.  The  multiple-hop  echo  structures  are  quite  similar  in  detail  to 
the  theoretical  ones.  Two  differences  can  be  noticed:  the  ordinary  and 
extraordinary  traces  can  be  distinguished  on  the  experimental  ionograms, 
and  the  noses  of  the  one-hop  echoes  appear  to  be  sharper  than  the  noses 
of  the  synthetic  ionograms.  These  differences  would  indicate  that  the 
actual  electron-density  profiles  were  not  as  blunt-nosed  in  the  vicinity 
of  the  electron-density  peak  as  were  the  mathematical  models  of  iono¬ 
spheres  used  for  synthesis. 

Figure  61  shows  experimental  evidence  of  the  presence  of  E  layers 
such  as  were  simulated  by  IID  012  and  089.  A  wide  variety  of  effects 
is  seen;  the  relative  strengths  of  the  E  and  F,,  layers  can  be  seen  to 
vary  from  one  record  to  the  next.  Figure  61b  shows  clearly  the  so-called 
"line  of  impenetrability,"  which  is  Holler's  term  for  the  locus  of  the 
left-hand  ends  of  the  F-layer  traces  in  the  presence  of  an  E  layer. 

This  line  can  clearly  be  followed  on  the  various  synthetic  ionograms, 
where  the  channel  between  the  E  and  F  traces  seems  to  follow  a  smooth 
locus  from  hop  to  hop. 
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Figure  61e  shows  striking  changes  in  the  relative  proportion  of  E 
and  F  echoes  at  Boulder  during  a  twilight  period. 

Figure  62  shows  five  records  that  are  similar  in  nature  to  synthetic 
ionograms  calculated  for  I  ID  133.  The  curious  shape  of  the  Fg  trace, 
which  comes  down  at  a  steep  angle,  turns  around,  and  goes  back  up  in  the 
direction  from  which  it  came,  can  be  understood  by  reference  to  the  133 
synthetic  ionograms.  These  patterns  apparently  reflect  an  ionospheric 
structure  that  is  closely  related  to  012  and,  in  fact,  the  transition 
between  the  types  is  evident  in  the  sequence  given  on  Fig.  61e. 


83  - 


SEL- 64- 106 


Best 

Available 

Copy 


Hi  T 

J  t! 

lull  i 

mil 

fiiiH  *  * 

IB  ii 

ijillrn 

piPM  1  •' 

gilt 

Ail  n 

!  |  '2|  1 

n  1 

X  H 

•  1  Ain 

tin  v 

>IU  m  ummm»  h 

ii  riu 

** 

B-rltr  i 

11 '  ii  iliU  i  * 

liAl 

rr. 

TlilSlB 

Hi  rib  ai 

1.1  ii 

in 

iiiiiii 

m  ii  if  ir  a 

:  I’ll 

IT 

i  mm 

i  l  - SI  lilt  Si 

11 

Ii 

mil  111  Hi 

««i  m  im 

,  111 

:n 

11  iilBB 

2 

■ 

1 

§ 

i 

i 

Bi 

Hi 

a 

2 

a 

2 

2 

2 

1 

iii 

8 

j 

u 

2 

2 

2 

S 

i 

B 

21 

ill 

8 

i 

m 

11 

2 

2 

2 

■ 

B 

Sj 

m 

1 

•2 

i 

2 

IS 

2 

■ 

a 

1 

■1 

ill 

u 

■ 

2 

2 

i 

8 

■ 

2] 

lag 

ii 

a 

1 

i 

■ 

a 

■ 

a 

21 

iii 

if 

l 

m 

FREQUENCY  (Me) 


From  Any,  Davies,  ami  Salaman  Kei  . 
record,  Sept.  12,  1952,  Section  I-l 
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OF  IMPENETRABILITY 


FREQUENCY  (Me) 


b.  From  Mo  Her  lief.  <1 


on  pane 


FREQUENCY  (Me) 


From  Any,  Davies,  and  Salaman  lb  1 
record,  April  2,  1952,  Section  I-l 
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FREQUENCY  (Mel 


FREQUENCY  (Me) 


1120  record 


d.  From  Agy  mid  Davies 
on  page  171. 


FREQUENCY  (Me) 


FREQUENCY  (Me) 


1 1  -1  <>  record 


From  Agy  and  Davies  Kel .  5 
on  page  I 


.  anil  N.i  I  aman  Hi  I 
I  952.  Sec  i l on  I -  I  . 


EXPERIMENTAL  lOXOGKAMS  SIMILAR  TO  THOSE  OF  I  ID  122 


e  io 

FREQUENCY  (Me) 


VIII.  SUMMARY 


A  method  has  been  described  that  permits  the  synthesis  of  artificial 
oblique  ionogrami:  by  means  of  a  modern  high-speed  digital  computer  at  a 
reasonable  cost.  Other  organizations  may  wish  to  carry  out  this  pro¬ 
cedure,  either  to  check  theories  concerning  the  origin  of  certain  oblique- 
ionogram  details,  or  to  build  up  a  more  comprehensive  "library"  of  well- 
understood  ionogram  structures.  Actual  computer  programs  are  not  in¬ 
cluded  here  because  the  potential  user  will  doubtless  prefer  to  receive 
such  programs  in  the  form  of  punched  cards. 

Thirty-nine  oblique  ionograms  are  synthesized  by  the  method  described 
in  order  to  illustrate  the  nature  of  the  ionograms  that  would  be  obtained 
experimentally  in  the  presence  of  ionospheres  similar  to  those  assumed. 
Simple  models  illustrating  the  most  gross  ionospheric  features  were  used 
for  this  first  effort.  However,  it  should  be  understood  that  the  method 
is  useful  for  any  electron-density  distribution  that  can  be  described  to 
the  computer.  This  includes  horizontal  gradients  of  a  wide  variety  of 
descriptions,  and  also  traveling,  localized,  ionospheric  irregularities. 

POSTSCRIPT  (December  1964) 

Prior  to  the  final  printing  of  this  report,  the  author  received  a 
private  communication  describing  similar  work  by  J.  W.  Finney  of  NBS  and 
H.  G.  Mol ler  who  was  temporarily  at  NBS,  They  made  only  one  ionogram, 
but  in  it  they  showed  the  effect  of  magnetoionie  splitting  for  a  partic¬ 
ular  path.  Readers  will  find  interesting  comparisons  between  our  tech¬ 
niques.  Used  together,  they  would  complement  each  other  quite  well, 
since  the  programs  given  here  are  fast  and  inexpensive  and  therefore 
suitable  for  extensive  investigations  or  preliminary  attempts.  The  more 
precise  (but  more  expensive)  method  of  Finney  and  Moller  could  be  used 
for  study  of  the  ionogram  perturbations  which  are  caused  by  the  geo¬ 
magnetic  field. 
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1  Attn:  Dr.  David  Johnson 

Mr.  Thurston  B.  Soisson 
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